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Propagation of light in a clumpy universe is examined. As an inhomogeneous matter 
distribution, we take a spherical void surrounded by a dust shell, where the "lost mass" in 
the void is compensated by the shell. We study how the angular-diameter distance behaves 
when such a structure exists. The angular-diameter distance is calculated by integrating the 
Raychaudhuri equation including the shear. An explicit expression for the junction condition 
for the massive thin shell is calculated. We apply these results to a dust shell embedded in a 
Friedmann universe and determine how the distance-redshift relation is modified compared 
with that in the purely Friedmann universe. We also study the distribution of distances in a 
universe filled with voids. We show that the void-filled universe gives a larger distance than 
the FRW universe by ~ 5% at z ~ 1 if the size of the void is ~ 5% of the Horizon radius. 

§1. Introduction 

One of the most interesting findings of the nearby redshift surveys is the discovery 
of large voids and that they appear to be a common feature of galaxy distributions. 
Great interest in such voids was first aroused by the discovery of the Bootes void,§) 
which was confirmed to have a radius of ~ 60 /i _1 Mpc.o At that stage, however, it 
could not be concluded that voids are common structures in the universe. The idea 
that voids are in fact common features of the large-scale structure of our universe was 
established by the recent surveys beginning from the CfA redshift survey,™ which 
introduced a picture of the universe where the galaxies are located on the surfaces 
of bubble-like structures with diameters 25 — 50 Mpc. Using the survey data such 
as CfA2§ and SSRS2,§ El- Ad et al.,% who established a method to find voids in 
these galaxy distributions, revealed that a substantial fraction of the volume of our 
universe is occupied by under density regions; they have shown that about half of the 
volume is filled with voids and that the voids have a diameter of at least 40/i _1 Mpc 
with an average underdensity of —0.9. 
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With the motivation of explaining such large voids, the gravitational evolution 
of a less-dense region has been studied by several authors. U'U> It was showncP that 
a less-dense region expands faster than the outer region, and a dense thin shell is 
formed behind the shock front by the "snow-plow" mechanism. The general relativis- 
ts motion of the void's shell was also studied bvJVlaeda and Sato.EP They adopted 
the metric junction method developed by Israel t^P and derived an expansion law of 
a void in the Priedmann-Robertson- Walker (FRW) universe. In addition to these, 
non-gravitational scenarios for the formation of voids have also been proposed. For 
example, in the explosion model, HIP a void surrounded by a thin shell is formed by 
the explosion of a pregalactic object. Another example is given by an inflationary 
model with a first-order phase transition, where vacuum bubbles are nucleated dur- 
ing inflation. These bubbles could represent the origin of voids. US' Thus the void 
model not only provides a simple picture of nonlinear density fluctuations but also 
is supported by many observational and theoretical considerations. 

There are also studies which investigate light propagation properties in a void 
system, such as the observational effect of voids on CMB anisotropics and the mod- 
ification of the redshift and luminosity for distant sources due to the existence of 
voids. The effects on the CMB anisotropics were discussed by several authors, H^'HiP 
suggesting that further CMB observations will give some constraints on the config- 
uration and origin of voids. Investigation of the light propagation properties in a 
clumpy universe using a spherical void model was undertaken by Sato,E3) who dis- 
cussed the modification of the redshift and luminosity in comparison with the FRW 
model. He found that the modification is third order in (Hr v /c) for the redshift and 
first order in {Hr v /c) for the luminosity, where H is the Hubble parameter, r v is 
the size of the void and c is the velocity of light. We believe, however, his result 
is suspect because the condition for the distance junction he adopted seems to be 
inappropriate. 

In this paper, we study light propagation in a void system based on the Ray- 
chaudhuri equation. A light ray bundle traveling across a spherical void surrounded 
by a thin shell suffers (de-)magnification and shear focusing. We calculate the modi- 
fication of the distance and redshift in comparison with the FRW model. Light prop- 
agation in the Swiss Cheese model has been studied extensively, t$' E_ZP but studies 
on a void surrounded by an overdense shell are far from complete. This paper gives 
an exact treatment of such a system. 

The organization of this paper is as follows. In §2 we derive the junction condi- 
tion for the expansion and shear of a null ray bundle across a shell. In §3 we calculate 
the modification of the distance and the redshift caused by the void-shell system and 
discuss its cosmological implications. The last section is devoted to a summary. We 
follow the signature of the metric and the convention of the Riemann tensor used in 
Ref. 18). 

§2. Basic equations 

In this paper, we consider a spherical void in the FRW universe. We assume 
that the interior of the void is a vacuum, so that the spacetime inside the void is 
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flat. The mass that would have been in the void is assumed to be distributed in a 
spherical shell surrounding the void. We treat the matter in the shell as an infinitely 
thin dust shell. 

In this section, we derive an expression for the angular-diameter distance in a 
flat spacetime and in FRW spacetime, which respectively correspond to the inside 
and outside spacetime of the void. We also derive the junction condition for the 
expansion and shear of a light ray bundle across the dust shell. 

2.1. Basic relations in geometric optics 

We use the notation used by Sasaki. @) Let k a be the tangent of a null geodesic 
parameterized by an affine parameter, A. We introduce a null vector field, l a , defined 
by 

l a l a = 0, l a k a = -l, l a . b k b = 0, (2-1) 

where the semi-colon denotes the covariant derivative. Using these vectors, we de- 
fine the metric H ab on the two-dimensional screen placed orthogonal to the spatial 
direction of the null geodesies: 

Hob = 9ab + kJb + lah- (2-2) 
The expansion 9 and the shear a of the cross section of the beam are defined by 

= k a . a , a ab = H a c {k c , d - k d , c - ^hJ] H b d . (2-3) 

We also introduce a dyad, {e A } (A = 1, 2), on the screen such that 

H ab = 5 AB e a A e b B , g ab e a A e b B = H ab e a A e b B = 5 AB and e% b k b = 0; (2-4) 

i.e. they are parallel-propagated alongthe geodesies. Then the basic equations of 
the geometric optics can be written asEI 

|J = -1244 - a AB a AB - \Q\ (2-5) 

ft B = -CA4B4 - 0O-AB, (2-6) 
dX 

with 

o-ab = o- ab e a A e b B , R^ = R ab k a k b , C A4B ± = C abcd e a A k b e c B k d , (2-7) 

where R ab and C abc d are the Ricci and Weyl tensor, respectively. Finally, we define 
the angular-diameter distance d\ by 

Tx^ = \°- < 2 ' 8 > 

By integrating Eqs. (2-5), (2-6) and (2-8), we can obtain the angular-diameter dis- 
tance. We can also write the equation for d\ as 

^d A jd K = ~R 44 - X -a AB a AB . (2-9) 
This expression is especially useful when the shear term vanishes. 
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2.2. Inside the void — flat spacetime 

In flat spacetime, the Raychaudhuri equations reduce to 

d „ 1 



n2 „AB 



and 



AX- r " (2 - 10) 



Suppose that the expansion 9 is positive initially, which implies that the derivative 
of c?a is positive. Then the angular diameter-distance is given by 



, A ^ AO { 1 + i, o( A-Ao)}|-{«^}, (2,2) 

where IcjaboI — \°~aboo~q B \ 1 / 2 ■ Note that when the shear vanishes, the distance is 
proportional to the time coordinate: 

d 

— dx = constant. (2- 13) 

By inspecting Eq. (2-12), one can see that the shear term appears as a correction to 
the shear-free case. 

2.3. Jumps of the optical scalars across the singular hypersurface 

The shell which surrounds the void generates a timelike hypersurface through its 
motion. Since we treat the shell as infinitely thin, the hypersurface may be regarded 
as singular.. We follow the procedure to treat a singular hypersurface formulated 
by Israel. HIP In order to obtain the jumps of the expansion and the shear of a null 
geodesic bundle across the shell, we proceed as follows. 

We consider a spacetime V with a timelike singular hypersurface S. The singular 
hypersurface £ divides V into two regions, V+ and V-. Consider a field variable W 
defined on E. The values of \P evaluated on either side of S will be denoted by )P± 
and the jump is defined by 

(2-14) 

We introduce a spacelike unit vector field n a (n a n a = +1) normal to U and the 
projection operator h% = 6% — n a nb- The extrinsic curvature K a b is then given by 

Kab = ~\h c a h d b £ n h cd = -h c a htV c n d , (2-15) 

where £ n and V a are the Lie derivative along n a and the covariant derivative, re- 
spectively. We have two extrinsic curvatures, K a b\±, which may differ, since U is 
singular. The jump of K ab across £ can be written in the form 

[K ab ] = 8ttG (s ab - \h ab si\ . (2-16) 



How Do Nonlinear Voids Affect Light Propagation ? 



5 



The tensor S a b here is related to the stress-energy tensor T ab through the equation 

S ab = lini v ( +£ T cd h c a h d b dx, (2-17) 

where x is a Gaussian normal coordinate (x = on £). Hence Sab is regarded as 
the surface stress-energy tensor of E. 

The jumps of 6 and gab on £ are calculated by integrating Eqs. (2-5) and (2-6) 

as 

r+e 

[9] = - lim / R 4A dx/k n , (2-18) 
[cjab] = - lim f +e Caaba dx/k n , (2-19) 

where 

fc n = k a n a . (2-20) 

Here we have set n a to point in the direction in which the affine parameter increases. 
From the above equations, we see that the jumps of 9 and gab on £ come from the 
diverging parts of -R44 and Caaba on U. Thus our next task is to calculate -R44 and 
Caaba- After a straightforward manipulation, we find 

i? 44 = (k n f£ n K c c + k a k h h c a h d h £ n K cd + finite term, (2-21) 

Caaba = { [e A e B - ~<Ub) k a k b - fc n (e^ + e\e\) k b + (k n ) 2 e a A e b B } 

x h c a h d b £ n K cd - ^S AB (k n ) 2 £ n K c c + finite term, (2-22) 



where e\ = e A n a . Using Eqs. ( |2T6 ) and (2T8)-(2-22), we finally obtain the jumps 
of 9 and gab as 

[0} = -^-k a k b S ab , (2-23) 



k n 

[<?ab] = ~- { (e A e B - Uab) k a k b - k n (e A e% + e B e\) k b + (k n ) 2 e\e b B 



x (s ah - h ab S^j - 27rG5 A Bk n S c c . (2-24) 

In this paper, we assume that the singular hypersurface which approximates the 
wall of the void has a surface stress-energy tensor of the dust type (we call it a dust 
shell), 

S ab = Du a u b , (2-25) 

where D is the surface density and u a is the 4- velocity of the shell {u a u a = — 1). In 
the case of the dust shell, Eqs. (2-23) and (2-24) become 

[9] = -^D(k») 2 , (2-26) 
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[(tab] 



4ttG 



-D 



2(e\k n - e\k n )(e B k u - e u B k n 



+ 



S AB {(k n ) 2 ~ (k u ?} 



(2-27) 
(2-28) 



where 

k u = k a u a and e\ = e a A u a . 

2.4. Outside the void — FRW spacetime 

We note here the analytical expression for the distance-redshift relation in flat 
FRW spacetime. The basic assumption is that the shear is negligible. We find the 
expression for d\ in terms of the cosmic time t as 



t x 2/3 
dA = dAi { 3 ( - 



ti 



2 - + ^d A i9i { 3 - - 3 - 



2ki 



2/3 * 



, (2-29) 



with an appropriate boundary condition for 8[, where k\ is the frequency of the ray 
observed by a comoving observer at t = ti. 



§3. Application and discussion 



3.1. Single void 

We consider a spherical void surrounded by a dust shell in a flat FRW universe 
(which we refer to in the following as a "FRW model"). The matter that would 
have been inside the void is distributed in the shell and we assume that the external 
universe is undisturbed. This is exactly the same as the model analyzed by Thomp- 
son and Vishniac.El) We consider a light ray bundle traveling across the void and 
calculate the redshift and the angular-diameter distance using the junction condi- 
tions obtained in the previous section. The result will be compared with the FRW 
relation. 

The metric of the FRW model is given by 

ds 2 = -dt 2 + a(tf{dr 2 + r 2 dQ 2 ). (3-1) 

The void is an empty region which can be described by a flat metric, 

ds 2 = -dt' 2 + dx' 2 + x' 2 dn 2 , (3-2) 

where the primes denote internal coordinates. The coordinate radius of the void is 
r v {t) in the external coordinates and x' v (t') in the internal coordinates. They must 
satisfy 

x' v (t) = a(t)r v (t). (3-3) 

Maeda and Satoi) and Bertschinger EIP showed that the shell radius expands asymp- 
totically as 

r v (t) tx tP with (5 ~ 0.13. (3-4) 
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In the following, the subscripts 1 and 2, respectively, denote quantities at the time the 
photon enters the void and at the time it leaves. Primes denote quantities measured 
by an observer at rest in the internal coordinate system, and unprimed quantities 
are measured by an observer comoving in the universe just outside the shell of the 
void. 
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Fig. 1. Cross section of a void and the photon path through the (j> = n/2 plane. The expansion of 
the void is not shown explicitly. 

The light ray bundle enters the void at t = t\ (corresponding to t' = t\ in the 
internal coordinates). Then the light arrives at the other side of the void at t = ti 
(t' = t^ 1 )- The arrival time is obtained by solving the equations 

t' 2 — t 1 = x'i cos 3q + x' 2 cos ^2 \ (3"5) 
x[ sin t^i' = x' 2 sin <^ 2 '- (3-6) 

(See Fig. 1 for the definition of the angular variables.) The redshift and the distance 
are obtained by solving the geodesic equation and the Raychaudhuri equations, as 
we do in the following. We also calculate the distance and redshift of a null ray in 
the FRW model which has traveled from t\ to t 2 . We then compare the distances 
and redshifts for these two models (see Fig. 2). 

Just before entering the shell, the light ray bundle has fci, d~ = &k\i #1 and 
o~abi = 0. At the time it crosses the shell, we impose Eq. (2-26) for 9, Eq. (2-27) 
for <7, and the following conditions for k and d^- 

k u a = k u a , k n a = k n a , cT^ = dj^. (3- 7) 

Note that these conditions for k a give the same result asthat obtained by the Lorentz 
transformation adopted by Thompson and Vishniacli^ From Eq. (2-12), the dis- 
tance at the time the ray bundle leaves the void is given by 

d A (t' 2 ) = dAi {l + ^r(t' 2 - t[)) x {1 + A}, (3-8) 
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Fig. 2. Trajectory of a dust shell and the path of a light ray in the void-shell system (left). For 
comparison, we consider a light ray which travels in the FRW universe for the same time interval 
measured outside the void (right). 



where A is the correction due to the shear term. This is to be compared with 
the distance in the FRW model given by Eq. (2-29). In the following analysis, we 
assume that the void is much smaller than the Hubble radius, and we define the 
small parameter rj by 

r] = a(t)r v (t)/t, (3-9) 

i.e., the ratio of the radius of the void to the Hubble radius. We expand all the 
quantities in powers of rj. 

Now let us consider the jumps of 9 and a. The jump of the expansion given by 
Eq. (2-26) is calculated as 

^ = gi-^/fci71 (1 + t,l00fl f l)2 ^ (3d0) 
1 9 J ' wi + cos^i ix' v ' 

where we have defined the velocity and the 7-factor as 

vi = a(h)^ = p Vl , 7l = (l- t ,2)-| j (3-11) 

and the ratio of the shell mass to the "extracted" mass as 

^Da 2 rl 3H 2 

* = (A li\ 3 3 Wlth ^ FRW = 'i - ^' ( 3 ' 12 

which will be set to unity, but at this stage we keep it as a free parameter for later 
discussion. We also used 



u a k a = 



— fci7i(l + v\ cos <Pi), n c k c = k\^i(v\ + cos^i). 



(3-13) 
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The shear induced by the shell is calculated using Eq.(2-27) as 

2 ^ l7l (l-t> 1 2 )sin 2 ^ 1 7 ?1 

= ^ = ~9 f , 1+ cos^ 7? (3 ' 14) 

where we have taken e' 2a = (0,0,0,3^). Since the increase in the affine parameter 
during the travel through the void is given by 

A - A = ^± = 2^- cos* m + 0( V j), (3-15) 

the correction term A in Eq. (3-8) is 

Ml „ab„ n x ,2 32 j2 (l-^)sin 4 ^icos 2 ^i 4 

A \~(J l cjabi(A-Ao =^rf ; ; -T-rz Vi- (3-16) 

81 {v\ + cos W\ Y 

Thus the correction due to the shear is of order rjf . 

Now we study the modification of the redshift due to the existence of the void. 
The redshift is defined in the void model by 

l + zy = ^ (3-17) 

and in the FRW model by 
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1 + zf=(^) 8 . (3-18) 



By a straightforward calculation, we find 
1 + z Y 



1 + Z F 



1 = (-^ cos 3 !£i + ^/Jcostfi) r,l (3-19) 



which agrees with the result of Thompson and Vishniac. Thus, the leading order 
of the modification of redshift is order nf. 

Next we study the modification of the angular-diameter distance. Expanding 
Eq. (3-8), we find 

<W(*2) = d A i + d A ip-0i cos "^t?! - d A i^-r]f 
k\ 9 

+d A i^#i (fl+l + l cos 2 $[) ii + 0( V f). (3-20) 



h 

For the FRW model, we find, from Eq. (2-29), 

dAF(h) = d>Ai + dAi r^di cosS^r/i - \ cos 2 

+d A ip-9t (P+I + I cos 2 &[) vl + 0( V f). (3-21) 
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FRW 



shell - void 



Fig. 3. Schematic picture for the volume of the matter which causes the Ricci focusing on the 
passing photon beam. In the FRW model (left), the beam is influenced by the matter contained 
inside the tube. In the void-shell sysytem (right), the shell contains only the matter inside the 
cones. 

Thus the difference reads 



Let us focus on the leading term in the difference. When cosl^{ is close to unity, 
i.e. when the beam passes near the center, the beam is de-magnified compared 
with the beam in the FRW model, to give a longer distance (we assumed / = 1). 
Contrastingly, when cos is close to zero, i.e. when the impact parameter is large, 
the beam is magnified. One can also note that the leading term of (g?av — <^af) for 
cos = 1 vanishes if / = 3. Recalling that / denotes the mass fraction of the dust 
shell to the mass that would have been in the void, we may say that in the void-shell 
system the Ricci focusing effect is weaker than in the FRW model by a factor of |. 
This can be explained in the following way (see Fig. 3). In the FRW universe, the 
beam suffers Ricci focusing by the matter within the tube shaded in the figure. On 
the other hand, in the void-shell system, the area of the shell which causes the Ricci 
focusing contains only the matter within the cones whose tops are at the center of 
the void. This explains the factor of g. 

Our result disagrees with the result obtained by Sato,© who argued that the 
difference is of order r\\. The disagreement comes from the junction condition for 
the distance at the time that the bundle crosses the shell. Sato assumed that the 
derivative of the distance with respect to the proper time measured by an observer 
at the boundary is continuous across the shell. However, we cannot read off any 
clear physical meaning that this condition should represent. Our condition, on the 
other hand, is obtained by integrating the Raychaudhuri equation, whose physical 
meaning is clear. As shown in the Appendix, the same result is obtained by a different 
method. Thus we conclude that the leading term in the difference is of order rj\. In 
other words, when a light ray passes only one void during its travel from the source 




(3-22) 
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to the observer, the modification of the distance compared with the whole distance 
is of order rfi, which may be considered to be very small. However, this means that 
the increase in the distance during the time that tye crosses one void is modified by 
the order of rji, since the radius of the void itself is of the order of r/i. Moreover, the 
effect is cumulative when the beam passes through many voids successively. Thus, 
it is interesting to study a multi-void system. This is the subject of the next section. 

3.2. Light propagation in multi-void system 

In this section, we study the distance-redshift relation in a universe filled with 
voids. The simple model used hereL_consists of randomly distributed voids, whose 
sizes evolve according to Eq. ( |3-4| ). tJ We start from a point source at t = t\ in the 
flat FRW universe. The light propagates in the FRW universe for a time interval 
~ 0.05 x t\ and then enters the first void surrounded by a dust shell. When the 
ray bundle crosses the shell, we calculate the expansion and the null vector using 
the junction condition described in the previous sections. We perform the same 
calculation for the time that the light leaves the void. For simplicity, we assume 
that the light ray enters another shell as soon as it leaves one void (see Fig. 4). The 
entering angle (or impact parameter) is treated as a random variable. At each exiting 
moment, we can calculate the redshift zy and the angular-diameter distance c?av 
from the source. We continue until zy exceeds a desired value. We also calculate the 
angular-diameter distance d\p for a source at the same redshift in the FRW universe. 
By repeating this sequence of calculations, we accumulate a large number of data 
(typically ~ 10 4 ) for dAv/^AF- We omitted cases when the expansion 6 becomes 
negative. The fraction is well below 1% even for the cases where the void size is 
about 0.1-ffJ -1 , and thus it has no influence on the following discussion. We also 
neglect the shear effect in the subsequent calculations, since the effect is negligible, 
as we saw in the previous section. 




, — *- observer 



Fig. 4. Photon beam traveling through the universe filled with voids. We assume that the beam 
enters another void as soon as it leaves one void. 

Figures. 5 and 6 display the probability distribution for c^av / ^AF • In Fig. 5 
we plot the distribution for samples at redshift between 1.0 and 1.1, and in Fig. 6 
for samples at redshifts between 3.0 and 3.2. The sizes of the void are adjusted so 
that the void has evolved at the observed time to radii of (approximately) O.OIHq 1 , 
0.05-ffJ" 1 and O.IHq 1 , where Hq 1 is the Horizon radius at the observed time. We 

*-* Our approach is similar to the method used by Holz and Wald.0 We have just substituted 
a void-shell system for their particle-hole system. 
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can see that the distribution is peaked above unity by 2% (4%) when the void size 
is OMH- 1 and by 4% (10%) when it is O.OSii^ 1 for z ~ 1.0 (z ~ 3.0), although the 
peak is rather broad for a large void size. The dispersion becomes even larger when 
the sample redshift is increased. Also note that the shape of the distribution is not 
symmetric. 



~ i 1 1 r 



z = 1.0 




< 



> 
< 

o 



0.5 - 



- 




d Av/ d AF 



Fig. 5. Probability distribution P of dAv/d\F for samples at redshifts between 1.0 and 1.1 . The 
normalization of P is such that the fraction of the beams which has dxv /d\F in the range [x, x + 
dx] is given by P{x)dx. The radius of the void at present is 0.01// ' 1 (square), G.Q^H^ 1 (cross), 
and 0.1// _1 (triangle) in each model. 



When the void radius is Q.1H~ 1 , the distribution seems to have a wavy feature. 
In order to see how this occurs, let us study higher order terms that we neglected in 
the previous section. Expanding Eq. (3-8) up to 0(r]f), we find 



0(vl) of 



^AV(^) 

h do 



- COS^ <?"l + 

O 



14 



4/3 2 + — p- - cos^i 



(3-23) 
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Fig. 6. Probability distribution of c/av I d\F for samples at redshifts between 3.0 and 3.2 . The 
normalization of P is the same with the previous figure. The size of the void is labeled by the 
same symbols as the previous figure. 



For the FRW model, we find 
0( V f) of 

do 

h d 

Thus the difference reads 

dAv(t' 2 ) - d A F(t' 2 , 



-3 (0+9 cos'** 



40 



14 



do 



2? cos 13 !fi + I 4/r + —p 



/16 t 1 d' 32\ 3r 3 

(27 + s iw cos M 



COS^i 



§(/ - 3« + if/ 



(3-24) 



cos !Pif7?. (3-25) 



This factor takes a small negative value when cos^i is small, i.e. when the beams 
passes close to the void wall. It decreases (the absolute value increases) as the impact 
decreases, then at some point it starts to increase, and ends up with a slightly positive 
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value when cos^i = 1. This causes the wavy distributions, which can be seen more 
clearly for a larger size of the voids. 

We find the peak of the distribution becomes close to unity when the void size is 
small. We can crudely interpretate this behavior as follows. While a lgiht ray bundle 
passes one void, the distance is affected by the amount of ~ 0(rf). The number 
of the voids which the lirhgt ray passes through before it reaches the observer is 
dA/r v ~ (1 + z)/r]. Thus, the accumulation of the effects amounts to 

AdA/dAj?(z)~(l + z)r}. (3-26) 

This consideration seems reasonable judging from the figures. 

Let us compare our model with the Swiss Cheese (SC) model. Lightnropagation 
in the SC model is studied using a similar approach by Holz and Wald. c3 According 
to their work, when the matter is condensated to point masses or compact objects, a 
large fraction of the beams from distant sources travel without being affected by the 
matter for most of their travel time. That is, a large fraction of the beams avoid the 
Ricci focusing effect, while a quite small fraction of the beams suffers strong shear 
focusing. Thus, the obtained distance for a fixed redshift favors a lower density 
model, even if the background model is a flat FRW model. On the other hand, in 
the void-shell system, a beam passes through a substantial fraction of the matter, 
resulting in a small difference from the FRW distance. In addition, when the beam 
successively passes through the voids, it is affected by almost the same amount of 
matter in the FRW universe as we have seen above. Thus, if the underdensity of 
the void is compensated by a surrounding (spherical) shell, the difference in distance 
from the FRW model is small. Also note that the distributions of distance are quite 
different for the SC model and the void model, reflecting the difference in the matter 
distribution. In principle, we can obtain information concerning the inhomogeneities 
from the scatter or distribution of the distance-redshift relation, though the current 
accuracy of observation is not sufficient to allow this. 

We finally note the implications of our results on the determination of the cos- 
mological parameters through the distance-redshift relation. As we have seen, the 
void structure increases the distance of a source compared with the distance in the 
FRW model; in this sense, the cosmological constant or a low density model is fa- 
vored. However, the magnitude of the modification is about 2% for a void whose 
size is 1% of the horizon radius, which is well below the current error in the estimate 
of the absolute magnitude (or distance) of the standard candles, such as a type la 
super nova. The scatter around the peak is also well below the current error as long 
as the void is not too large (< 1% of the horizon radius). 

§4. Summary 

We have studied light propagation properties in a spherical void (empty region) 
embedded in a spatially flat FRW universe. The void is surrounded by an overdense 
wall, which we treated as a dust shell. We considered the null geodesic equation and 
the Raychaudhuri equation, and derived the junction condition for the expansion and 
the shear across the shell. Using these, the redshift and the distance were calculated 
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in the void-shell system, and we compared the results with the FRW relations. We 
have found that the modification of redshift by a single void is of order tj 3 , and the 
modification of the distance is of order rj 2 , where r\ is the ratio of the void size to 
the Hubble radius. 

We have also discussed the cosmological implications of voids by considering a 
universe filled with voids. The void- filled universe slightly increases the mean dis- 
tance and produces a dispersion in the distribution of the distance. This may work 
to lower the density parameter or favor the cosmological constant when the cos- 
mological parameters are estimated through the distance-redshift relation of distant 
sources. However, for samples with redshifts ~ 1.0, the magnitude of the effect is 
well below the current error in estimating the distance of the standard candles as 
long as the voids are within the observationally favored size. 

Appendix A 

Simple Derivation of d\ for a Single Void 



Our analytic method based on the geometric optics equations (2-5) and (2-6) is 
not only useful for the present analysis but also can be extended to general space- 
times including singular hyper surf aces. For the model of a void in the flat FRW 
background, however, we have a simpler derivation of the angular-diameter distance 
which we now present. 

Thompson and Vishniac derived not only the redshift deviation (3.19) but also 



the scattering angle of a photon: 



PL3D 



5a = -<Zq + + &2 - $2 = -vi sin(2^ 2 )- (A-l) 

9 

As depicted in Fig. 7, we define a as the angle formed between the direction of 
observation and the direction of the void's center, Dl as the comoving distance of 
the void's center, and D$ as the comoving distance of the light source. Hereafter the 
subscript S denotes a quantity at a light source. In the two-dimensional comoving 
space we denote each position by the vector symbol r and introduce a vector basis, 
[n(a),l(a)]. The trajectory of a photon between a source and an observer was 
obtained in Ref. 14). It is given by 

*s ~ ro = 3t (1 - Vas)n + (l - ^/a^ - -^j Sal 

= D s n + D LS Sal, (A2) 
where Dl$ = D$ — Dl. In the absence of a void, this relation reduces to 

(rs - ro)b = 3*o (1 - \/as)n = D s n, (A-3) 

where the subscript b denotes a background unperturbed quantity. 

The angular-diameter distance of a light source is originally defined as 

d A = j, (A-4) 




Fig. 7. Cross section of a void on the ip = tt/2 plane of the model. We depict the trajectory of a 
photon from a light source to an observer. Define a as the angle formed between the direction of 
observation and the direction of the void's center, 5a as the scattering angle of a photon, Dl as 
the comoving distance of the void's center, and Ds as the comoving distance of the light source. 
We denote each position by a vector symbol r, and introduce a vector basis, {11(a), 1(a)}, in the 
two-dimensional comoving space. 




where T> is the proper distance across the source and 5 is its angular diameter. Here 
we calculate the manner in which D is changed by a void for a fixed 5, considering the 
radial direction from the void's center and the circumferential direction separately 
(see Fig. 8). The radial component T>^ and the circumferential component T>^ are 
written as 



V (r) 



d(rs - r Q ) 



da 



Aa, 



p( c ) oc y = (r s - t ) ■ e y , 



(A-5) 
(A-6) 
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where y, e y and Aa are denned in Fig. 8. A straightforward calculation leads to 

1 " D S da ' (A ?j 

1 = J" S 5a cot a. (A-8) 

Thus, the order of magnitude of the modification of the distance is essentially deter- 
mined by the scattering angle 5a, which is of order r/|- 
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